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Abstract. We present an explicit cosmological model where inflation and dark 
energy both could arise from the dynamics of the same scalar field. We present our 
discussion in the framework where the inflaton field <f> attains a nearly constant velocity 



in 



\d(p/dN\ = a + f3exp(f3N) (where N = In a is the e-folding time) during inflation. 
We show that the model with |a| < 0.25 and (3 < can easily satisfy inflationary 
constraints, including the spectral index of scalar fluctuations (n s = 0.96 ± 0.013), 
tensor-to-scalar ratio (r < 0.28) and also the bound imposed on during the 
nucleosynthesis epoch (fi^lMeV) < 0.1). In our construction, the scalar field potential 
always scales proportionally to the square of the Hubble expansion rate. One may 
thereby account for the two vastly different energy scales associated with the Hubble 
parameters at early and late epochs. The inflaton energy could also produce an 
observationally significant effective dark energy at a late epoch without violating local 
gravity tests. 



PACS numbers: 98.80.Cq, 98.80.-k, 95.36,+x [arXO706.2654 
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1. Introduction 

The WMAP measurements of fine details of the power spectrum of cosmic microwave 
background (CMB) anisotropics [1] have lent a strong support to the idea that the 
universe underwent an inflationary expansion in the distant past [2] . The WMAP data, 
along with the independent observations of the dimming of type la supernovae in distant 
galaxies [3] also favour a result of growing evidence that a large fraction of the energy 
density of the present universe is 'dark' and has a negative pressure, thereby leading to 
the ongoing accelerated expansion of the universe. It is then natural to ask whether it is 
possible to unify the inflation and quintessential fields. In a viable theory the primordial 
inflation may lead to have a dark energy effect in the conditions of concurrent universe. 
This picture merits broader discussion. 

The main observation that has led many to believe that the dark energy is Einstein's 
cosmological constant A, for which w\ = p A / p A = — 1 identically and at all times, is the 
concordance of different cosmological data sets, which appear to indicate that the dark 
energy equation of state w DE = P DE /p DE is not much different from —1 at the present 
epoch. This solution to dark energy, however, raises two immediate questions: (i) why 
is pa = A/8ttG ~ 3pM today? and (ii) why is p A (~ 1CT 12 (eV) 4 ) so tiny? Apparently, 
p A is fifteen orders of magnitude smaller than the electroweak scale (m EW ~ 10 12 eV), 
the energy domain of major elementary particles in standard model physics, and it is 
not known at present how to derive it from other small constants in particle physics. 

The cosmological constant as the source of dark energy is only a possibility. The 
other possibility is that the cosmological constant (or gravitational vacuum energy) is 
fundamentally variable. Explicit examples are provided by models that use a dynamical 
scalar field with a suitably defined scalar potential V(<f>). Quintessence models 
are among the most popular alternatives to Einstein's cosmological constant as they 
generally predict at late times a small (but still an appreciable) deviation from the 
central prediction of Einstein's cosmological constant, i.e. w A = —1. Observations only 
require that w DE < —0.82 at present epoch [HE], so one finds worth studying models 
that support a time-varying dark energy. 

There are arguments in the literature [U [5] that an appropriate modification of 
Einstein's theory provides an alternative resolution to dark energy problem and a natural 
framework to address the inflationary paradigm. In this context, higher-dimensional 
braneworlds models, scalar-tensor theories and R + f(R) gravity models, which derive 
motivations from the original idea of Kaluza and Klein to its modern manifestation in 
string theory, have been of particular interest. 

A simple modification of Einstein's theory of general relativity, which involves a 
fundamental scalar field <fi with a self-interacting potential V((p), is given by 



where k is the inverse Planck mass nip] = (StiGn) 1 ^ 2 and C m is the matter Lagrangian. 
This theory has been studied over the last three decades by crafting different types of 
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scalar potentials. The list of the potentials can be frustratingly long, which includes 
the quadratic potential V{4>) = \ m 2 ^ <p 2 widely considered in inflationary contexts and 
the inverse power-law potential V(4>) oc <p~ a (a > 2). These examples are perhaps 
sufficiently simple to understand the basic ideas of inflation and/or the dynamics dark 
energy in the concurrent universe, for a review, see [6], but they hardly explain the 
cosmic expansion of our universe exhibiting all relevant cosmological properties. It is 
thus natural to ask whether it is possible to unify the inflation and quintessential fields 
by finding (or constructing) a more general potential. 

The model of quintessential inflation [7] proposed by Peebles and Vilenkin uses 
the idea that inflaton potential could end up as an effective present-day cosmological 
constant [8] or quintessence [9]. Although quite appealing, the potential considered 
in [7j, which consists of two parts: V(<f)) = A(0 4 + M 4 ) (0 < 0) for inflation and 
V(<f)) = AM 8 / (0 4 + M 4 ) (0 > 0) for quintessence, finds no natural field theoretic 
motivations. Recently, attempts have been made in constructing a working model of 
quintessential inflation within the context of higher dimensional braneworld models, see, 
e.g. P21 [11] and references therein for the earlier proposals. Also, there are suggestions 
that a unification of the inflationary era (triggered by R 2 type corrections) and the 
late-time acceleration can be made through a simple construction of the modified F(R) 
models [5], as well as within the framework of reconstruction of scalar-tensor gravity [I]. 

In this paper, we reconstruct an explicit observationally viable model for evolution 
from inflation to the present epoch by maintaining the structure of the theory defined 
by ([I]). Our reconstruction approach yields a smooth, exponential potential that 
describes both the inflation and quintessential parts. The model can be shown to be 
compatible with current cosmological observations, and, presumably, it can be embedded 
in higher dimensional theories of gravity, such as string theory. 

The rest of the paper is organized as follows. In section 2, we motivate our 
construction with an appropriate ansatz for an inflaton field. We then invert the 
system of autonomous equations to determine the inflaton potential, along with other 
cosmological variables. There we also find conditions that have to be satisfied by 
the reconstructed potential to be consistent with the WMAP inflationary data. In 
section 3, we briefly discuss about an efficient method of reheating, so called 'instant 
preheating', applicable to our model. In section 4, we include the effect of ordinary 
fields and then find an explicit quintessence potential in a background dominated by 
radiation (or matter). In section 5, we show how the reconstructed potential produces 
an observationally significant effective dark energy and its associated late-time cosmic 
acceleration. In section 6, we discuss on a possible way of evading local gravity 
constraints imposed on the model. Further generalization of our construction with 
higher-order corrections is briefly discussed in section 7. Finally, section 8 is devoted to 
conclusion. 
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2. How might inflaton roll? 

In this section, we neglect the effect of ordinary fields (matter and radiation). The set 
of autonomous equations of motion following from (JT]), with C m = 0, is given by 

2" 



V{4>) = m 2 H 2 



3 - 2m 2 



1 alH 



(2) 



2 fldH\ 

where H = H(<f)) = a/ a is the Hubble expansion parameter and a(t) is the scale factor 
of a spatially flat Friedmann- Robert son- Walker universe. 

One of the most crucial parts of a consistent inflationary model is to understand 
the time-evolution of the inflaton field 0. Any choice of should give rise to a flat 
potential as required for inflation and also be consistent with cosmological observations, 
including WMAP results. To this aim, a simple (and possibly a natural) choice for the 
evolution of inflaton field is 

= O - a m p ln[a/ aj - 7 m p (a/ 2C , (4) 

where |a| < 0(1) and a i is the initial value of the scale factor before inflation. We shall 
take 7 = 1 for a reason to be explained below, while the parameters a and £ (< 0) will 
be fixed using bounds on inflationary variables inferred by the WMAP observations [I]. 
The evolution of the inflaton field in (jlj), or equivalently 0(t) = 0o + colnt + c\jt v 
(with p > 0), is a generic solution for a modulus and/or dilaton field in many four- 
dimensional string models, see, e.g. [12]. The assumption (jlj) holds, almost universally, 
in many well motivated inflationary models that satisfy slow roll conditions, after a 
few e-folds of expansion. For instance, for the chaotic model of inflation with the 
potential V(0) oc m 2 2 , one has a oc e^ 2//2 (cf equation (2.4), ref. [13]) and thus 
|0| = v^21na + const. As discussed in |14| . even for two scalar fields model, if the 
slow-roll conditions 3if0^ ~ V<p. are satisfied at Hubble exit, then M = In a depends 
linearly only on the field values, leading to a generic situation that 0(t) oc lna+ (small 
correction) . 

The reconstructed scalar field potential is given by 



V(0) = m^ 2 (0) (3 - % (0)) , (5) 



where 



H{(j>) =Mexp 



a 2 C 
-^iV(0)-aX-|x 2 



fdH/d(j)\ 2 1 2 



(0) = 2m^— =-(a + 2(Xy, (6) 

where X = 7e 2 ^ N ^', N((f>) = lna(0(t)) + C. Note that the parameter 7 appears 
only in a combination with e 2 ^; so using a shift symmetry in and/or choosing 
the constant C appropriately, we can always set 7 to unity, thus 7 = 1 henceforth. 
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The energy scale M (which appears as an integration constant) can be fixed by 
the amplitude of density perturbations observed at the COBE experiments, namely 
(dV/d(l))- 1 V^ 2 /(y/7Eirml ] ) ~ 1.92 x 1(T 5 . With a = 0.2 and W e = \n(a f / ai ) ~ 55, 
assuming that ( < 0, we find M ~ 7.4 x 10~ 5 m p = 3.1 x 10 14 GeV. 

With a slowly varying e H ((p), the scalar curvature perturbation can be shown to 
be [15] 

P}l\k) = T-^-^-(l - e H r^ ( , (7) 



r(3/ 2 r *> \2^\ /aH __ k 

where ^ = 3/2 + 1/ (p — 1) and a oct p . The scalar spectral index n s of the cosmo logical 
perturbation is defined by 

B .<*> sl+ -SEif- (8) 

The fluctuation power spectrum is in general a function of wave number k and 
evaluated when a given comoving mode crosses outside the horizon during inflation: 
k = aH = a e H((p)e~^ N is, by definition, a scale matching condition. Instead of 
specifying the fluctuation amplitude directly as a function of k, it is convenient to 
specify it clS cL function of the number of e-folds. 

In the case a = 0, we get H(4>) oc exp (/>{2(j) — 4>) and e H (<f>) = 2( 2 K 2 ((j) — (f) ) 2 . 
The scalar potential takes a familiar form: V((j>) oc m 2 [3-2CV(0-0 o ) 2 ], where 
m 2 oc -f/" 2 . The number of e-folds is M c = -J= (e^)" 1 / 2 ^ = ^ In , where 
4> 2 < <fi l < 4> Q . Since r] H = ^ (d 2 H / d<f> 2 ) / H = e H — 2( is small only for a limited 
range of inflaton values, <fi ~ 0o ; the number of e-folds is large only when ( is very 
small. In this case, however, almost no gravitational waves would be produced, leading 
to an exponentially suppressed (close to zero) tensor-to-scalar ratio. The spectrum of 
scalar (density) perturbations is also almost Harrison- Zeldovich type, n s = 1. This last 
result is, however, not consistent with WMAP observations [1]. Thus, without loss of 
generality, we demand that \a\ > 0; more precisely, 

C < 0, - 2(e 2CN <a<VQ 

so that V{4>) > 0. The spectral index n s is now given by 

, n A o? + 6a 2 \ + 12aA 2 + 8(A + 8A 3 
n -l = 2r)„-&e„ = , (9) 

(up to leading order in slow roll parameters) where A = (e 2 ^ c and Af e is the number 
of e-folds of inflation between the epoch when the horizon scale modes left the horizon 
and the end of inflation. 

The scalar spectrum on scales accessible to CMB observations is perhaps that 
measured at the instance when observable scales exit the horizon during inflation. In 
most models this corresponds to a phase of inflation between e-folds 50 and 60. We 
summarize the results in a Table (for M e = 50 and £ = —0.1): 
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Figure 1. The tensor-to-scalar ratio r ~ vs the scalar spectral index n s , with 

a = 0.21,0.20 and 0.19 (top to bottom) and C = (-0.2,0). The solid (dotted) lines are 
for We = 60 (W e = 40). 
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Figure 2. The scalar spectral index n s vs a, with Af e = 70 (solid line), J\f e = 40 
(dotted line) and C = {~ 2,0}. Except for a < |C| < 0.05, the value of n s does not 
much depend on 





n 
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r = l6e H 


a 




r < 0.28 


n s > 0.965 




< 0.18 


< 0.017 


n s = 0.96 




0.32 


0.200 


0.020 


r = 0.1 


0.987 




0.112 


0.006 



In figured] we show the plot of tensor-to-scalar ratio r with respect to n s , and in 
figure [2] the plot of n s with respect to a. Within our model, both n s and r do not much 
depend on the number of e-folds except when ( is positive, which we reject on physical 
grounds. 

With the WMAP3 bound on the tensor-to-scalar ratio, r < 0.28, we find n > 0.965 
f° r C ^ —0.1. The bound r < 0.28 implies that e H < 0.0175 and imposes a 
relation (for a given N) between A and a. Using equation (Q, we get n s > 0.965 
for ( < —0.1. The scalar spectrum is red-tilted except in the case that a < and 
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Figure 3. The running of scalar spectral index, a , with respect to a and £ with 
J\f e = 60. Except for a < |£| < 0.05, a s does not much depend on the number of 
e-folds. 



both ( and r are sufficiently close to zero, e.g., for ( = —0.005 and r = 0.001, we get 
(a, n s , N) = (-0.0051, 1.0107, 50), (-0.0057, 1.0097, 60). 
The running of spectral index, a s , is given by 
dn s dn s dN d<p 
as = dink = dN d<p dink' ^ 

where 

_d^_ _ y/2^(0) dN 1 

d\nk~ mp {l-e H {<P)Y mp d<P~ y/tojfi' 1 J 

These relations hold independent of our ansatz (j4j). In our model, the value of ct s is 
found to be small, when satisfying 0.01 < a < yz and ( < (cf figure [3]). 

We conclude this section with a couple of remarks. Studies in [16] show that, 
in slow-roll inflation, one may relate the variation of the inflaton in terms of e-folds 
N = ln(a//aj) to the tensor-to-scalar ratio r 

— ^ = — = xF- (12) 
mp dN m P y 8 

The WMAP bound on the tensor-to-scalar ratio is r < 0.28 (95% confidence level). This 
then implies that a < 0.187 in the present construction. This is completely consistent 
with our discussion above. 

The reconstructed potential may be expressed as 

V (fP) = K^T [ 6 " (« " 2 CV? " 2«C , (13) 



where 



H(if) oc exp 



a 2 C 
-jN(<p)+a<p-i(<p + aN(<p)) r ' 



(14) 
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Figure 4. The scalar potential for some representative values of a = 0.3, 0.4, 0.5 
(top to bottom), C = - 0.1 and C = -10. 



where N(ip) = ln[a(tp(t))} and tp = (</> — o )/mp. The shape of the potential (as 
depicted in figure Hj) as well as its functional form is qualitatively similar to a class of 
scalar potentials one would obtain via warped flux compactifications of string theory, 
see, e.g. [TTj . The predicted characteristics of inflationary phase (of the potential) can 
easily be made to comply with the WMAP results pQ. So our method of reconstruction 
may be considered as a point in favour of providing a believable physical basis for the 
inflation. Moreover, a large part of our construction does not depend on the details 
of string theory or the dynamics of scalar fields abundant in any higher dimensional 
theories but has a general validity, and thus would remain useful even if string theory 
is invalidated. 



3. Reheating after inflation 

A satisfactory model of inflation should perhaps be followed by a successful reheating. To 
this end, the 'instant preheating' mechanism presented in [18] and applied to exponential 
potentials in [19] might perhaps be the most efficient method for reheating the universe. 
Here we briefly outline a viable mechanism of reheating in our model, leaving the details 
for future publication. 

According to (jlj), after a few e- folds of inflation, since (N < 0, one has <ft ~ 
—a m p H. Clearly, with a < \^2, the kinetic term never dominates the potential term. 
As a result there remains the possibility that the expansion enters inflation from which 
it never recovers. So our model has a chance to work only if the matter and/or radiation 
energy density terms sometime after inflation is large enough to dominate the inflaton 
energy density. 

Without loss of generality, we can make the inflation end at the origin by translating 
the field 

V(<p) = M 2 m 2 p (3 - a 2 jl) e a ^/ m p) + small correction, (15) 
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so after inflation ip = (0 — en d) < 0. Following [TSJ [TH] we assume that the inflaton 
field p interacts with another scalar field x- The interaction Lagrangian is 

Ant = -t^W - hipipx, (16) 

where g and h are coupling constants, and ip is a Fermi field. The production of x 
particles commences when the adiabatic condition 

\m x \ < m\ (17) 

is violated, i.e. when \p>\ > gp 2 , where m x = g\tp(t)\. So, the particle production may 
occur when 



, I < Vend Q 1/2 Knd 4 (g) _ ,„ / 1R N 

l^l^y^ gX/4^1/2 = *w ( 18 ) 

The process of particle production occurs nearly instantaneously, within the time 

prod i^J ( } 

during which the field <p remains in the vicinity of tp — 0. As the field rolls to <p < 
direction, the mass of the x particles begins to grow, since m x = g\tp(t)\. The occupation 
number of x particles is n k ~ e _7r(kA Vod) 2 

, with k being; the canonical momemtum. The 
energy density of particles of the x field created in this process is 

Px = rn x nJ°^)\ (20) 



where the number density n x = (2-it 3 )^ 1 J °° k 2 n k dk ~ (27r)~ 3 (ag) 3 / 2 V^(p). If the x 
particles can rapidly decay into fermions or the quanta of the x field were to convert (or 
thermalize) into radiation, then the radiation energy density would increase sharply to 

5/2 3/2T/3/4 

Pr*fbc~ g7r 3 ^ V pr od ~ 0-0027/ a V en ». (21) 
At the end of instant preheating 
P v 370 

Although p x / p v is small quantity to begin with (for any generic value of the coupling 
g < 0.3 and a < V^), P x ( or the decay product of the x field) will decrease as a~ 3( - 1+w ^ 
(w < 1/3) and come to dominate p^ since the field tp is rolling down an exponential 
potential and its energy density could decrease much faster p oc 1/a 6 after inflation. To 
illustrate this one considers a cosmic evolution by suppressing dV/dp, so <p + 3Hp = 0, 
whose solution is <p = <p + p l j a~ 3 dt. According to (jlj), p ~ —am p H and hence 
a(i) = a end (t + (3/a)t) 1/3 . We then find 

a 2 „, 4 _„_ Q , /a -\ e 



2 ^y^)~lO^ 9 m 4 p ^j . (23) 

For a < \/3) there would be no kinetic regime. Nevertheless, since p x (or the decay 
products of the x field) may decrease much slower l/[a(t)] n (n < 4) than p^, it will 
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eventually dominate the scalar energy density before the production of light elements or 
the BBN epoch. Instant preheating may be followed by reheating which occurs through 
the decay of x particles to fermions as is evident from the interaction term in ( fTBT) . 



4. Growing matter 

Given that the inflaton field decays to some radiation and heavy particles, it would 
be natural to expect, at later stages of inflation, small but nonzero values for both 
the matter and radiation energy densities. The growth in matter energy density can 
naturally affect (or modify) the form (or shape) of the scalar potential, leading to an 
additional term in the potential with a relatively large slope. This last feature is perhaps 
required to make our model compatible with the big-bang nucleosynthesis (BBN) bound 
imposed on the scalar field energy density. 

Here we take the matter Lagrangian in its simplest form, which is Einsteinian 

C m = C(g^i ^m) = V^9 (P M + Pr) , (24) 

where p,^ oc a~ 3 ^ 1+W( ^\ i = M (matter) or R (radiation). Of course, one could allow 
in principle an explicit coupling between the 0-field and matter. It is believed that 
inflation was followed by an instant preheating (or reheating) and then by a radiation 
dominated phase, so the strength of coupling between the field and matter could be 
neglected during both the inflationary and the radiation-dominated epochs. Any such 
couplings, however, can be relevant at later stages of evolution, especially, at galactic 
distance scales (see section [6]). 

The set of autonomous equations of motion that follows from equations (OQ) and 
( 124|) may be given by [20] 



K 2 V{(f>) 



(3 + e)(l-fl tt ) + in; 



2 ~ w 

..2 a' 2 



# 2 (0), (25) 



=Q^-2e(l-Q w ), (26) 

= -2ea w -3(i + w)n w , (27) 

where Q w = Qm + Qr, the prime denotes a derivative with respect to N = ln[a(t)], and 

During radiation dominance would remain small but nonzero. This last assumption 
is consistent with the fact that the fixed point solution Q w = 1 is always unstable. 
Notice that the behaviour V(<f)) oc H 2 (<f)) holds also in the presence of ordinary fields 
(matter and radiation). 

From equations (|25l)-(l271). along with equation (j4j), we find 

«V(0) = [3(1 - n w ) - 5(0)] , (29) 

e(0) = + 5(0), (30) 
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" C + 3(1 +w)J(-C( ( j ) ))dN( ( j ) y lJ J 

where C is an integration constant and 

B(d>) = \ (a + 2Ce 2 ^)) 2 , (32) 

C(0) = e ( Q ^ 3(1+w) ) NW exp [2a e 2 ^ + Ce 4CNW ] . (33) 

As compared to the inflationary potential given by equations (1131) and (114)) . we now 
have the effect of matter fields (matter and radiation together). Of course, in the limit 
that Q w — > 0, equation (T29T) reduces to ([TBI . 

During radiation domination, since Q R S> Vt M and Q, ~ Q R , we have w ~ 1/3. 
One also notes that the last term in equation (j3J) does not contribute (significantly) after 
inflation. Therefore, from equations ( [29l) - ([3~Tl) . we get 

^ = 1 + w - « 2 /3 ffg e 3(i +w ) re (^-^)/ a ^ , 34 ^ 
™ 1 + w H 2 (4>) 

where H 2 (<j)) = #2[ e <*«(^i) + e 3(*+i)K(*-fc)/a] ^ and 

V(0) = S [ aie 3(i+w) K (^-0 2 )/a + a2 e«-W-^)] , (35) 

where ct\ = «2 = ^f-j an d -^o, 01, 02 are integration constants; we take 

< 2 <C Exponential potentials of a such form, which also arise ubiquitously in 
particle physics and string theory models [2T] , by themselves are a promising ingredient 
for building a natural model of quintessential inflation. In order for the scalar field 
potential not to dominate the energy density of the universe during BBN, it is required 
that 3(1 + w) > y/Ea, which easily satisfies the bound imposed on during the 
nucleosynthesis epoch, fL (1 MeV) < 0.05 [22]. By taking a < 0.8 and w ~ 1/3, we 
correctly reproduce a double exponential potential anticipated in [23] . 

The reason why the quintessential part of the potential, equation (135]) . has a 
different form with respect to its inflationary part is easy to understand in our model. 
During inflation the matter contribution (and its possible coupling with the inflaton 
field) can be safely neglected. This is, however, essentially not the case for quintessence 
part. Another source of this difference is that the last term in equation (j4j) does not 
contribute (significantly) at later stages of evolution, like during the radiation-dominated 
epoch. 

5. Late time acceleration: Dominance of dark energy 

At late times, without loss of generality, one takes p R p M and w ~ 0. One 
also assumes that is rolling only slowly, such that \<j>\/H < m p . In this case the 
inflaton potential takes a simpler form, for the evolution of the universe could naturally 
lead the potential part to dominate the kinetic part: 2V(0) oc 2 , with m being the 
proportionality constant. Explicitly, we get 

y(0)=m 2 p i/(0) 2 -^ T (l-fiJ, (36) 
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Figure 5. Evolution of the universe passing from matter dominance (w c s ~ 0) to 
scalar field dominance (u> e ff < —1/3), with m = 3,5, 10 and 50 (from top to bottom). 



e((j)) = — ^-j- — | wQ w and Q w = 1/[1 + S(z + l)~ 3ui ], where z is the redshift factor and 
w = w + The Hubble parameter H(cj)(z)) (and hence V(^(^))) can be expressed 
in a closed form using the relation e = H/H 2 . The numerical constant 5 can also be 
fixed using observational input: an ideal situation would be that the universe re-enters 
into an accelerating phase (e > —1) for z < 1. The universe passes from a decelerating 
phase to an accelerating phase when Q w < ™~\ ■ The dark energy equation of state is 
w 4> = P(f>/ ' P<t> = i+2m ' therefore, with m = 50, we get w c s = —1 — 2e/3 ~ —0.76 and 
w, ~ —0.98 (see also figure [5]). 

The behaviour of dark energy similar to that depicted in figure [5] may be seen 
directly from equations ()34j) - (j35p . Using the relation e N = e lna = (1 + and making 
the assumption that ordinary matter (including cold dark matter) is approximated by 
a non-relativistic perfect fluid and p R <^ p M , so that w ~ Pm/Pm ~ ^' we 

n ~ n M = f i - — 1 - 1 — (37) 

V 3 y i + c (i + 2)« 2 - 3 v ; 



and 



e(^) = -1 - q{z) = -\n M - ^, (31 



where q is the deceleration parameter. Hence 



H(z) = H 



n m0 (l + ^) 3 + c (l-fi m0 )(l + ^r ' • (39) 



1/2 



The numerical coefficient c Q may be fixed such that Q M = 0.27 at z = 0. With a < y/2, 
the second term on right-hand side decreases with z at a slower rate as compared to 
p M (which varies as (1 + z) 3 ) as well as to that of the curvature, p k (which varies 
as (1 + z) 2 ), so naturally exhibits 'dark energy' as late times. As depicted in 
figure E], for a ~ 0, the universe starts to accelerate when z < 0.8. For a larger a, 
acceleration starts at a lower redshift; with a moderate value of a ~ 0.26, we get 
w DE = w 6 ^ (a 2 -3Q 6 )/3Q^ -0.97. 




Figure 6. The deceleration parameter q(z) with respect to redshift z, and a = 
0.8,0.6,0.3,0.01 (top to bottom). The free parameter c Q in equation ([37|) is chosen 
such that Q M0 — 0.27. 

The present model addresses the cosmic coincidence problem, only partially. In 
fact, the cosmic coincidence problem (i.e. why p^ ~ 3p M now?) often involves some 
kind of fine tuning, and it is not an exception here. An interesting observation is that 
this last phenomenon requires either a specific ratio between the kinetic and potentials 
terms, or a specific value for the field velocity = K(p/H, which is characterized by the 
parameter a, so as to realize a quintessence dominance for z < 0.85. 

6. Evading gravity constraints 

In the above discussions we ignored the coupling of the 0-field with matter. This is 
perhaps justified. 

The dark energy or the cosmic acceleration problem is essentially a problem 
associated with largest cosmological scales: in order for the field to play a role of 
dark energy its effective mass should be at least in the range of the present value of 
the Hubble parameter, H Q ~ 1CT 33 eV. In turn, one takes the runaway quintessence 
potential satisfying ^/V((j)) ~ H ~ 1CT 33 eV; the range of the interactions mediated by 
the scalar field can be of the order of the Hubble horizon size. However, Newtonian 
tests of Einstein's general relativity and fifth force experiments such as the Cassini 
satellite experiment put stringent bounds on the gravitational coupling of light scalar 
particles. That means, a putative dark energy field should be sufficiently massive at 
much smaller scales. Thus a mechanism similar to that in Chameleon field theory [24|, 
which combines both a quintessence-like behaviour leading to dark energy at late time 
and a gravitational coupling to matter which is appreciable in high density regions, 
could be operative in our model. 

To this reason, one allows a nontrivial coupling between the 0-field and matter, 
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and, accordingly, takes the matter Lagrangian in a general form 

C m = Cty m , A 2 (0)<^) = v^A 4 (0)^p w , (40) 

where p^ oc a~ 3( - 1+Wi \ a = aA((f>). couples to the trace of the matter stress tensor, 
g^tTfij , so the radiation term p R does not contribute to the equation of motion for 

f>4> + 3H P<j> i 1 + W d) = -fa<xM)PM> ( 41 ) 



where V = (1 - 3<), = ±0 2 + V ((/>), = pjp^ w. = pj Pi and a<j> = 
Equation (14ip . along with the equation of motion for ordinary fluids 



p\ + 3^(1 + Wi) = +0r/a A(0)p;, (z = M,R), (42) 
guarantees the conservation of total energy, namely p + 3H(p + p) = 0, where 
P = Pa/ + Pr + /V 

In the discussion below we take w M = 0. The effective scalar potential is then given 

by 



V eS = V{4>) + p M J a^A((j)) dcj), (43) 

where p M oc 1/a 3 . For \aj > 0, the model needs to be confronted with the present-day 
equivalence principle bound, a 2 < 5 x 1CT 5 . On largest scales probed by WMAP, where 
p M ~ p crit ~ 10 -12 (eV) 4 (where p crit = 3H 2 /8ttG n is the critical energy density), the 
last term above is only sub-leading, which is suppressed by a factor of a^. In turn, 

can be sufficiently light, m, = V^ 2 ~ 10~ 33 eV ~ (10 28 cm) _1 , and its energy density 
may evolve slowly over cosmological time-scales. But within solar system distances, 
where p M is roughly 10 23 times larger than its value on large (Hubble) scales, the term 
proportional to p M can be more relevant. On Earth, p M ~ 10 30 x p crit , the Compton 
wavelength of the field can be sufficiently small, A c ~ m" 1 ~ 0.1 mm as to satisfy local 
tests of gravity. That is, in high density (and high curvature) regions the quintessence 
field may end up almost in a squeezed state. 



7. Further generalization 

Although the model above is canonical in describing the basic ideas of quintessence, 
there exist theoretical and phenomenological motivations for studying modifications of 
the Einstein-Hilbert action which allow non-trivial couplings of to some quadratic 
Reimann invariants (of the Gauss-Bonnet form 1Z 2 = R flu \ a R IJ ' uXcr — AR^vR^ + R 2 ) and 
antisymmetric tensor fields [25], [26] 

C = v ^{~& + " ~ Q ^ n2 ) + (44) 

where £(0) = - ±(<90) 2 -V(0), H 2 = H v ^W vX and H vX = d [fM B wX ] is the antisymmetric 
3- form field strength. Allowing Q{<p) ^ in f|4*4l) . one introduces a pseudoscalar degree 
of freedom a, via the ansatz H^x = ^fg e^ v \ T d T a. Like 0, the axion field a is a function 
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only of time. In particular, the coupling J 7 ((f)) allows new cosmological solutions for 
which the dark energy equation of state can be less than —1. To be precise, we note 
that 



k (Pde+Pde) 



K 2 (f)' 2 + (1 - e)£V + (45) 



where fi^ = 8JFif = QF'H 2 . The antisymmetric 3-form field does not modify this 
equation because it contributes to p DE and p DE with the same magnitude but with 
opposite signs, namely, K 2 p DE /H 2 = x 2 /2 + y 2 + SQ^ + 3Q g and n 2 p DE /H 2 = 
x 2 /2-y 2 -(2+e)^-^-3fi g , where tt g = 2Q(<f)) a' 2 , x = K<p/H and y = kVV/H. We 
can get = p^/p^ < — 1, without requiring a superluminal expansion e = H / H 2 > 0, 
or having to introduce a non-canonical (phantom) field. Most features of the model 
(CQ) would arise in the limit where J-~((f))7l 2 and Q((f))7i 2 are sub-leadings to V((f>) (see 
below). 

In the above model, the axion field a does not play much role. With a generic choice 
of (?($) = £ e 2 * (where $ = (f)/m p ), the 5-field equation of motion, V M (e 2$ #^ A ) = 0, 
is solved for H^ uX = e~ 2 *e M1/Ar 9 T cr. The integrability condition, d^H^xr] — 0, yields 
a + 3#cr + 2<M = 0. With the ansatz (g]), we get 

^ oc exp [(2a - 3)N + 2(e 2 < N ] . (46) 

After a few e-folds of inflation, the last term above would become small, since (N < 0. 
The scalar potential reads 

Vffl = ^-[6 -a 2 - 12Q <WW] , (47) 

where H = H exp[a<&/2]. This result in conjunction with equations (j4j) and (146]) implies 
that for a < 2, Q((f))Ti 2 decreases faster than the scalar potential V((f>). 

Next we briefly discuss some qualitative features of the reconstructed scalar 
potential with a nonzero T '((f)). With the ansatz (j3J), and with G((f>) = 0, the 
reconstructed potential is given by equation fl5J); the parameter £ H ((f>) reads 

1 



e H ((f)) = -(a + 2(Xy + 3n :F 

= \(a + 2(X) 2 - 24H 2 (a + 2(X)^^-. (48) 
2 d(p 

Clearly, in the case < 0, the coupling J 7 ((f)) could increase the period of inflation by 

making e H smaller. This effect can be opposite in the case > 0: it could be that 

inflation ended due to a slowly increasing derivative of the coupling, dJ 7 /d(f>, such that 

~ 1/3. 

With J 7 ((f)) 7^ 0, the corresponding potential may be reconstructed by providing 
an extra condition or by demanding a specific relation between the functions V((f>) and 
J 7 ((f)) (see, e.g. [27], where a general method of reconstruction was developed, including 
the effect of scalar-Gauss-Bonnet coupling). In the particular case that a(t) ~ a e H o*, 
we find 

n -- e —- a 2 -^e 2 ^-^e^ (49) 
^ 3tf 2 1-2C 1-4C ' 1 J 
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where N((j>) = lna(0(t))+ const. Again, after a few e-folds of inflation, since 
exp[2CiV(»] < 1, we get 



e N(4>) V ^ 5 e *(0 

^ = -«-3^' # = 1+ 2 a+ 3^- (50) 

This result reveals a generic situation that the coupled Gauss-Bonnet term is only 
subleading to the potential V(<p). This behavior of our model may be present also when 
the Hubble parameter changes appreciably with e-folding time, as happens at later 
stages of inflation. 

The presence of ordinary fields (matter and radiation) in our model does not 
introduce much complication, apart from slightly modified expressions for V(4>) and 
for the added degrees of freedom come with additional equations of motion. 



8. Conclusion 

We have presented an explicit cosmological model for evolution from inflation to the 
present epoch that we believe satisfies the main observational constraints, including fine 
details of the power spectrum of cosmic microwave background anisotropies, e.g., a red- 
tilted scalar spectrum with small tensor-to-scalar ratio, r < 0.28, the bound imposed on 
during the nucleosynthesis epoch and present epoch local gravity tests. It is therefore 
potentially of great utility. 

In our analysis, just one assumption, equation (j3J), that is regarding the evolution 
of inflaton field, has been made, which is indeed a common feature of many motivated 
slow-roll type inflationary models. Moreover, for a slowly rolling inflaton field, m p <f)' = 
m pJw < 0.25, the gravity waves or the amplitude of tensor perturbations can be 
suppressed in our model. This might actually be needed in our model, in order to 
satisfy the BBN bound. 

The present proposal also simplifies the role of the inflaton by almost decoupling 
it from the (background) matter on large cosmological scales. On the scale of the 
solar system, due to the large surrounding matter density, the dark energy field can be 
sufficiently massive, e.g., ~ a/ A c fr, ^ > 10~ 3 eV, thereby quenching the deviations 
from Einstein's gravity on distances larger than a fraction of millimeter. Moreover, the 
model possesses an attractor behaviour for the inflaton and matter densities analogous 
to the tracking solution of, e.g., the inverse power- law potential, V(4>) oc <p~ a with a > 2. 

The model proposed here may provide a reasonable explanation to the question: 
why is the cosmological vacuum energy smalV. The interpretation of gravitational 
vacuum energy (or dark energy) in our framework is likely to yield V(4>) < 3(1 — 
Q mQ )H 2 m p and exhibit scaling behaviour for p^, being proportional to the square of 
the Hubble rate. As a result, within our model, p^ ~ 2H 2 m 2 p = 2 x 10~ 66 eV 2 m p ~ 
3.5 x 10 _47 (GeV) 4 would be the most probable value of dark energy density at the 
present epoch. 
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n s = 0.96 ± 0.01 and r ~ O(10~ 3 — 10~ 4 ), see also [28] for other details. This result is 
compatible with WMAP5 data [29] . 
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